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Abstract 

We set up a trace formula for the relativistic density of states in terms of a topolog- 
ical sum of classical periodic orbits. The result is applicable to any relativistic integrable 
system. 
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1 Introduction 

Gutzwiller's trace formula of 1971 expresses the density of states g{E) of a quantum me- 
chanical system approximately as a sum over all periodic classical orbits JI||. Later, Balian 
and Bloch || presented a formula which also applies to nonintegrable systems. It arose 
from a study of sound spectra in cavities with reflecting walls of arbitrary shape in two and 
more dimensions. Gutzwiller's formula applies only to systems with isolated orbits. It fails 
if there exist degenerate families of periodic orbits connected by continuous symmetries [[|. 
The problem arise in the derivation of Gutzwiller's formula from a stationary phase approxi- 
mation to the trace integral over the semiclassical Green function at fixed energy. It contains 
an oscillating exponential of the eikonal function S"(x, x; E) = § p • cfx of the periodic orbits 
passing through the point x. A continuous symmetry makes this independent of x over an 
entire spatial region swept out by the symmetry operations. Then the second derivatives of 
the eikonal function vanishes in that region, resulting in a divergence of the stationary-phase 
integral. Strutinsky and coworkers |4j] removed these divergences by going back the convolu- 
tion integral in the time- dependent propagator and performing exactly as many integrations 
in that integral and in the trace integral over the Green function, as there are independent 
parameters describing the degeneracy. Later, Creagh and Littlejohn || pursued the same 
idea in a generalized phase space which also contains room for the continuous symmetry of 
the system. For integrable systems, their procedure is similar to that of Berry and Tabor || 
who derived a trace formula for integrable systems employing the action-angle variables. 

All this development has so far been restricted to the nonrelativistic regime where the 
particle solves the Schrodinger equation in some external time-independent potential. The 
purpose of this paper is to begin adapting the methods to relativistic particles described by 
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the Klein-Gordon equation in the external potential. Our final result will be a relativistic 
generalization of Gutzwiller's trace formula, expressing the density of states as a topological 
sum over the relativistic closed classical orbits. The formula is applicable to integrable 
relativistic classical systems. 

Relativistic quantum mechanics is of course not really a consistent theory. At rel- 
ativistic velocities, particles will be created and absorbed, and the particle number is no 
longer conserved, thus violating the current conservation law of the Klein-Gordon equation. 
Quantum field theory is certainly the appropriate tool to describe relativistic particles. In 
the classical regime, however, the particle number is fixed and these problems are absent, 
so that a semiclassical expression for the density of states in terms of relativistic classical 
periodic orbits is a consistent approximation expected to render a reliable results for those 
systems in which particle creation and annihilation play only a minor role. 

2 Relativistic Quantum-Mechanical Trace Formula 

Consider a relativistic particle of mass m in an external time-independent potential V(x), 
whose quantum mechanics is governed by the Klein-Gordon equation 



where c and H are speed of light and Planck's constant, and A(x) is a magnetic vector 
potential. We have shifted the energy origin to the rest energy mc 2 in order to have a smooth 
limit to nonrelativistic bound-state energies. Since the potentials are time-independent, the 






Et / n ^(x), and (H) takes the Schrodinger-like 



form 




(2) 
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where 

E 2 — m 2 c 4 . . 

and He is the Hamilton operator 

H E = p 2 /2m + [2EV(x) - y 2 (x)]/2mc 2 , (4) 

with p = —ihd x . It is useful to view as a special case of a more general eigenvalue 
equation 

W E *(x) = f *(x), (5) 
which arises from a Schrodinger-like equation 

W B *(x,r) = ina r *(x,r) (6) 

by a factorized ansatz \&(x, r) = e~ l£T ^ h ^f{x). Then the variable r plays the role of a 
pseudotime, and the Hamilton operator He is the pseudotime-evolution operator governing 
the T-dependence of the system. 

Let ^ n ( x ) be the eigenf unctions of Eq. (|) with eigenvalues £e(n). Then the physical 
energies E n of the particle are given by those values of E at which the pseudoenergy is equal 
to e: 

£ En (n)=e. (7) 

As an example, consider the Coulomb potential V(r) = — e 2 /r of the relativistic 
hydrogen atom. Equation (E) leads to the radial eigenvalue equation 



d 2 R(r) 2dR(r 



+ r-^ + 



<ir 2 r cfr 



2m / Ee 2 /mc 2 \ 1(1 + 1) 
h 2 \ E r I r 2 



a 



R(r) = 0. 
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Its solutions yield the bound state pseudoenergies depending on principal quantum number 
n and angular momentum /, but degenerate in the azimuthal quantum number m: 



E 2 /mc 2 



£ E (n,l,m) 



a 



2 ■ 



n — 1; 2, 3, • • • , s 
i = 0, 1,2, ••• ' {) 



(n-l- 1/2) + + 1/2) 2 - a 2 
Inserting these into Eq. (|7|), we obtain the well-known relativistic bound energies of the 
Coulomb system: 



E n ,i = ±mc 



a 



-1/2 



(10) 



{n-l- 1/2) + yj(l + 1/2) 2 - a 2 
The complete information on the spectrum of eigenvalues of the Klein-Gordon 
equation ([I]) is contained in the pole terms of the trace of the resolvent R(E) = 
i(e - He{^) + it])^ 1 : 



r{E) =Tr R(E) = i Tr [e - H E (n) + ir]]-\ 



where the infinitesimal positive quantity rj guarantees the causality of the time dependence 
of the Fourier transform of r(E). The imaginary part of r(E) defines the density of states: 



g(E) = -Imr(E) = Ti5{e - H E ). 



(12) 



In terms of the eigenvalues £ E (n), the density ( P^) has the spectral representation 



g(E)=Y,5(s-£ E (n)), 



(13) 



where the sum over n covers all quantum numbers. This sum will now be performed in a 
semiclassical approximation as a sum over periodic classical orbits. 

For the sake of generality, we assume that the particle moves in D-dimensions, and 
assume that the motion has been transformed to D cyclic coordinates whose motion can 



easily be quantized (torus quantization) . The labels n will then be integer- valued vectors n = 
(ni,n 2 , • • • , rip) with non-negative components rij. For the purpose of deriving a semiclassical 
approximation to (jig) , we convert each sum over rij = 0, 1, 2, ... in Eq. ( |l3l) into an integral 
with the help of the Poisson summation formula [[7| ^] 

OO OO 

£/(») = £ / (i4) 

n=0 fc=-oo 1/0 

Here we have assumed that the function f(n) and its derivatives with respect to n vanish at 
infinity, and the lower limit 0~ on the integral sign indicates that the integration starts on 
the left-hand side of the origin to include the entire ^-function generated by the sum over k. 
The superscript will be omitted in the sequel. Thus we obtain 

g(E) = J2fd D n6 (e-S E (n)) e 2 ™ kn , (15) 
k J 

where each component of the integer-valued vector k = (ki, k 2 , ■ ■ ■ , ko) runs form minus to 
plus infinity, while the now continuous variables rij are integrated from 0~ to infinity. 

For integrable systems, the integration variables rij in Eq. (|15"1) can be replaced by 



the values of the action integrals appearing in the relativistic quantum conditions [10 



Ii = — (f p • <ix = (rii + ^ ) h, (16) 



2?r JCi V 4 

where p is the relativistic momentum of the point particle along closed loops Ci on an 
invariant torus. The quantum numbers are the same nonnegative integers as above, while 
/ij are the numbers of conjugate points along the orbit Cj. Thus we can rewrite Eq. (|T3]) as 

j /*oo poo POO 

9(E) = T7T£ e ~ ikW2 / dh dl 2 --- dI D 5(e-£ E (I))e 2M , (17) 

where we have changed the argument of £e(n) to £e(I), and introduced vectors /j, = 
(/ii, /i 2 , . . . , Hd)- 
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Consider now the lowest term with k = 0, for which the oscillating exponentials in 
Eq. (|17[) are absent. It contributes a smooth density of states 

j rco rco roo 

g(E) = - u dh dl 2 --- dI D 5(e-S E (I)), (18) 
ri Jo Jo Jo 

where the lower bounds of the integral has been moved to zero, since the classical orbits 
for k = have zero length, making /x equal to zero. The multiple integral (|T8|) is just the 
classical density of states 

UE) = / / d D pd D q5(e-£ E (p,c l )), (19) 

which in cyclic coordinates reads 

j roc r2n roc r2n roc rlix 

9d{E) = — — ^ / dh d(px dl 2 d V2 --- dI D d(p D 6(e-S E {I)), (20) 
{27rn) u Jo Jo Jo Jo Jo Jo 

reducing to (18) after integrating out the angular variables. The classical density of states 

is also referred to as the Thomas-Fermi density [[J, 

We now turn to the oscillating k ^ parts of g{E). With the help of the integral 

representation for the S-f unction 

5 (e-£ E (I)) = — - / dre^-^ 1 ^, (21) 

we rewrite this as 

S g (E) = -L r dr-^ Y> e -W«l* [°° dh r dl 2 --- H d/ De ^M-(i)P 

2irh J-oo h ^ Jhni/i Jhn 2 /4: iv D /4 

(22) 

where the primes on the summation symbols indicate the omission of k = 0. The integrals 
over h and r are now evaluated in the stationary phase approximation. Let us abbreviate 
the action in the exponent by 

A k (I, r) = 27rk • I+r[e-£ E (l)}. (23) 



Its extrema lie at some I = I, r = t, where 



9A 



t f - 9t 

1=1, T = T 



I=I,T=f 
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The first set of equations yields the semiclassical quantization condition 



(24) 



27rJfei = fWi(I), i = 1,2, (25) 

where 

<*(I) = (26) 



at I are the angular velocities for the pseudoenergy e. The solutions of Eq. (P5[) yield actions 
I as nonlinear functions of k and f: 

I = I(k,f). (27) 
From Eq. (^5] ) we obtain the important relation for the resonant tori 

£ = '•./ i. >.•••./;. (28) 

Since fcj ere integer numbers, the orbits on the torus must have commensurate frequencies, 
so that only closed periodic orbits contribute to the density of states in the saddle point 
approximation. This establishes the connection between 5g(E) and the relativistic periodic 
orbits of the classical system. If the frequencies are not commensurate, the orbits do not close 
although the motion is still confined to the torus. Such orbits are called multiply periodic 
or quasi-periodic. 

Each relativistic periodic orbit is specified by k; it closes after k± turns by 2tt of the 
angle ipx, k 2 turns of <f2, ■ ■ ■ ■ Thus k plays the role of an index vector characterizing the 
topology of the periodic orbits. For this reason, the sums in Eq. ( |T7D is also called topological 
sum. Note that Eq. (|28| ) admits only /^-values of the same sign. 



The second equation in (p3[) specifies r via 

e - S E (l(f (£))) = 0. (29) 

Having determined the saddle points, the semiclassical approximation requires the calcula- 
tion of the effect of the quadratic fluctuations around these. For this we expand Eq. (|23|) up 
to the quadratic terms, and shift the integration variables from I to I' = I — I. The lower 
bound of the integrals is then transformed into hfi/A — I. For sufficiently large actions Jj, the 
sharpness of the extrema at small % allows us to move the lower bounds to minus-infinity. 
This approximation is excellent for highly excited states. We now perform the Gaussian 
integrals and obtain the oscillating part of the relativistic density of states 

5gW(E) = l-^Jh D+1 ^' e -»W2 e -W4 _L_ |detM|f 1/2 e« hrk - i /' i , (30) 

^ k ^~ 

where M is the stability matrix 



M = 

V di, 

whose determinant is, according to formula 



/ _ d 2 £ E d£ E \ 

T 

dlidlj dli 



(31) 



det ( ^ B j = det A det(D - CF A~ X B) (32) 



given by 



det M = det H ^H^iv. (33) 



where 



H = (34) 

The Maslov index v is equal to N + — N~ — iV°, where iV ± denote the numbers of positive and 
negative eigenvalues of matrix Hij, and iV° is unity (zero) if the sign of u T if _1 u; is positive 



(negative). The second factor in ( |32"D has been simplified using the equation of motion for 
the cyclic variables (p: 

^ = Vi&(I) = «C0, (35) 



the right-hand side being also equal to 



2vrk , . 

u> = — . (36) 



Since for every k there is an equal contribution from — k, we may replace the exponential by 
a cosine and obtain 



5g {2) {E) = ^-J27i/h D+1 5T _ m , W9 det H uj t H' 1 uj _ 1/2 cos [2vrk • (i/ft - p/l) - tiv/A . 

27T — ' x' I L \ / 

k 

(37) 

The relativistic trace formula fl3~T|) gives us a basis for understanding quantum phenomena at 
the relativistic level in terms of classical orbits. In general, we just need to evaluate the clas- 
sical Se flj for integrable systems, and consider some shortest orbits. As in nonrelativistic 
systems, we expect astonishingly accurate energy spectra from Eq. (|37|). 



3 Three-Dimensional Relativistic Rectangular Billiard 

As a first application, consider the motion of a relativistic particle in a three-dimensional 
rectangular billiard with sides of length a lt a 2 , and a 3 along g 1; q 2 , and q 3 axes. The quantum 
spectrum of Eq. (^) with Dirichlet boundary condition is given by the pseudo energies 



£EK,n 2 ,n 3 ) = ^- + + ^ = 1,2,3) = 1,2,3, (3* 

ZTYI \ CL-^ Q>2 ^3 / 

The physical relativistic energy spectrum is obtained from Eq. (|7|): 



\ 



2 2 2 



7T 



2 ^ 2 (^ + 7f + 3l+^ 4 - (39) 

til (X 
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As in the nonrelativistic case, this result is exactly reproduced by the relativistic quantization 
according to Eq. (jlBD . The numbers //j are all equal to 4, since the wave functions have 
Dirichlet boundary condition. At every every encounter with the wall, the action picks up a 
phase it. The relativistic action variables are therefore 



h = — f Pidqi = Uih, 2 = 1,2,3; ra = 1,2,3, 

Z7T 



(40) 



The classical Hamiltonian may be expressed as 



s E (i) =— 4 + ^ + -| , 

2m V ag / 



(41) 



and the corresponding angular frequencies are 



7T 



(jj; 



ma 



i — 1)2,3. 



(42) 



We now determine the saddle points I. According to Eq. (^5|), these are given by 



(/i,J 2 ,/ 3 ) (r) 



2ma\ki 2ma\k<i Ima^k^ 



TIT 



TIT 



TIT 



(43) 



leading to the pseudoenergies at the saddle point 



2m 



£ b (I(t))=— EC** 



(44) 



The saddle-point value of r is determined by (p9|), yielding 



'2m 



r 



A 



i=l 



(45) 



From these saddle point values, we obtain 



2e, 



1(f) 



(46) 
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so that the sign of UiH i j 1 Uj is positive and the number N° in the Maslov index v = N + 
N~ — N° vanishes. The determinant of the second-derivative matrix is 



d 2 £ E _ vr_ 6 

dlidlj m 3 afa^ai, 



All eigenvalues of the matrix 8 2 Se/ 'dlidlj are positive. Thus we identify the indices N + = 
3, N~ = 0. Inserted into Eq. fl37|), we finally obtain for the oscillating part of the relativistic 
density of states 

where jo(x) is the spherical Bessel function of order zero jo(x) = sin(x)/x. The symbol L 
denotes some lenght scale which may be any average of the three length scales ai, 02, or (13, 
while 



denotes the energy associated with L. The quantity S (k) is 



S (k) = -VE 2 - m 2 c 4 2\jk\a\ + k\a\ + k\a\ = p L k . (50) 

It is precisely the relativistic eikonal pL^ of the classical periodic orbits of momentum p 
and total length L^. In general, the inclusion of only a few shortest orbits in Eq. ( f48D 
yields the correct positions of the quantum energy levels. The three-dimensional relativistic 
rectangular billiard may serve as a prototype of the relativistic semiclassical treatment for 
arbitrary billiard systems. 

Let us compare the calculation of (|48|) from our relativistic trace formula ( |37|) with a 



direct calculation from an inverse Laplace transformation of partition function Z(j3), i.e, 

1 re+ioo 

g(E) = — dPePZffl, (51) 



2-rri 



e—ico 
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where the partition function is given by 



oo oo oo 



Z(P)= E E E ex V {-p£ E ( ni ,n 2 ,n 3 )}, (52) 

ni=l n2=l ri3=l 

with the pseudoenergies Se^i, n 2 , n 3 ) of Eq. fl38|) . The problem is the same as in the 
calculation of the Casimir energy for the box. Since ( |52"D is a product of three independent 
sums 

oo 

W)= E ex P{-^(^)}> ^ = 1,2,3, (53) 

n;=l 

we may process each sum separately Applying the Poisson formula (O) to the sum over Ui 
we find 

Z0) = E / dme-^i^Me 2 ***" 1 * - - 

ki=—oo 00 — 



Inserting this into fl52p, and using the integral formula 



e^e-^= - JJ2V^ , (55) 



2m Je-ioo \K, 

we obtain the exact level density of the relativistic three-dimensional rectangular box 



5 *W. 

(56) 

The leading term comes from a proper three-fold sum, and is given by 

with S (k) of Eq. fl5"0|). This agrees with the semiclassical result (f|£|). The second set of 
terms gives corrections from the faces of the box: 

ki,k2=— oo \ / 
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where 

S 2 (k u k 2 ) = -^E 2 - m 2 c 4 2y/kfal + k\a) = pL klM (59) 

(2) 

are the eikonals of the orbits on the faces. The functions g\j (E) are the level densities of 
the planar facial "boxes" . 

The third set of terms in (|56| ) stems from the edges of the box, being the level density 
for these one-dimensional "boxes" of length af 

^-sfcj>»PF)- (60) 



where 



1 



Si (k) = -VE 2 - m 2 c 4 2ka % =pL klM . (61) 

These boundary terms can be obtained also from the general trace formula ([14]) by calculating 
higher-order corrections to the semiclassical approximation (|37|). 

The last term in (]56f) is a delta function at e = which does not contribute to the 
level density at e > 0. 

The classical (Thomas-Fermi) contribution to the density of states is 

9{E) = E~ [l)jE-Tz -82L* + 8 V ~ L ) " (62) 

Here V3 = 010203 is the volume of the box, V 2 =2 (a\a 2 + a 2 az + 0103) the total surface, and 
V\ = 01 + 02 + 03) the sum of the edge lengths. 

4 Concluding remark 

For relativistic integrable systems, we have derived a semiclassical trace formula by trans- 
forming the relativistic quantization conditions into the topological sum involving all closed 
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relativistic classical orbits. Certainly, our final result ( pTj) can also be obtained by an ab ini- 
tio procedure, starting out from the relativistic path integral for the relativistic fixed-energy 
amplitude representation ]7|, |TT], [12 



G(x b , x a ; E) = JjL j™ d L J Dp<$> [p] J D D xe 



DiA E /h 



(63) 



with the action 



rn 

Ae [x, i] = dr 

J Tn 



M 



2p(r 



e 
c 



-x 2 (r) + -A • x(r) + 



2Mc 2 



(64) 



T"6 



where L is defined by 

L= / ' drp(r), (65) 

with p(r) being an arbitrary dimensionless fluctuating scale variable, and is some con- 
venient gauge-fixing functional, such as <3> [p] = 5 [p — 1]. The prefactor Ti/Mc in ( |B"3|) is the 
Compton wave length of a particle of mass M, the field A(x) is the vector potential, V(x) 
the scalar potential, E the system energy, and x the spatial part of the D + 1 -dimensional 
vector x = (x, ir). This path integral forms the basis for studying relativistic potential 
problems. Choosing p(r) to be equal to unity, the amplitude (B3) becomes 



G(x b ,x a ;E) = j^J o dLexp -eL J V D xexp 



-A, 



(66) 



where the fixed-energy action Ae is given by 



A, 



io IT X ( T ) + - A - X ( r ) 



2Mc 2 



\/ 2 (x) - 2EV{x) 



(67) 



The semiclassical approximation to the relativistic fixed-energy amplitude ( |S6]) is |13 



G sc (x 6 ,x a ;£) = ————— — V / dLe ieL/h 



class. traj. 
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x det 



-d x id a jA E (xi,,Xa',L) 



1/2 



,iA B (x6,x a ;L) /%—i-KV /2 



(68) 



The associated density of states is obtained from the trace of (§H): 



class. traj. 



x / cPxdet 



-9 a .i9 a ,MB(x 6 ,x a ;L) 



1/2 



,iA E (xi,,x a ;L)/h—iTrv/2 



(69) 



The trace operation in Eq. (pq) is integration over all periodic orbits in the pseudotime "L" . 
If the relativistic systems is integrable, it can be expressed in terms of action-angle variables 



as 



d D x det 



-d xl d xi A E {y. b ^ a ;L) 

V r d D ipL- D ' 2 det 
^ Jo 



1/2 



,iA E {x.b,x- a ;L)/h-inu/2 



d 2 S E (I) 
dlidlj j 



1-1/2 



i[2irI-U.-S E (l)L]/h-iiTv/2 



(70) 
(71) 



A E (x 6 , x ; L) = I • ( <p b - <p a ) - £ E {l)L = 2ttI ■ k - S E {I)L, 

thus establishing contact with the earlier treatment in which r plays the role of L. 
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